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Abstract. Bauer and Itzykson showed that associated to each labeled map embedded on an oriented 
Riemann surface there was a group generated by a pair of permutations. From this result an algorithm may 
be constructed for enumerating labeled maps, and this construction is easily augmented to bin the numbers 
by the genus of the surface the map is embedded in. The results agree with the calculations of Harer and 
Zagier of 1-vertex maps; with those of Bessis, Itzykson, and Zuber of 4-valent maps; and with those of 
Ercolani, McLaughlin, and Pierce for 2i/-valent maps. 

We then modify this algorithm to one which counts unoriented maps or Mobius graphs. The results in 
this case agree with the calculation of Goulden and Jackson on 1-vertex unoriented maps. 



1. Introduction 

The recent work of Ercolani, McLaughlin, and Pierce 6 analyzed the fine structure of families of prob- 
ability measures on the space of N x N Hermitian matrices. The principle measures analyzed in [6] are of 
the form 

(1.1) dfH = exp (-NTv [Vt(M)}) dM, 

where dM is the product measure 



dM = ( J] dRe(M y )cflm(M ii ) dM h 

l<i<j<N I \l<i<N 



the potential is 



V t (X) = ^X 2 + t\ 2 



and where Zm is the normalization of the measure. The question of interest in that paper was the asymptotic 
structure of log(.Zjv) for large values of N, and its analytic dependence on t. 

More generally, the partition function of random matrices, introduced above as the normalization factor 
of the probability measure 1|1.1|) . is defined as 

(1.2) Z N (t 3 ,t 4 ,...,t 2u ) = J exp^-NTi(^M 2 +t 3 M 3 + t 4 M 4 + ---+t 2ly M 2, ^dM, 

where the integral is taken over the space of N x N Hermitian matrices. The partition function has been 
used extensively as a model for the partition function of 2D-Quantum Gravity |2|, EH- 

Derivatives of Zttits, i/t, ■ ■ ■ , tiv) evaluated at tj = capture the Gaussian moments J H Tr(AP) fc d/io(-M)- 
These moments are evaluated with the Wick lemma in terms of pair correlators of the matrix entries. In 
analogy to Feynman diagrams, there is a correspondence between labeled maps and the terms in the Wick 
expansion. This fact led to the conjecture that the logarithm of the partition function (|1.2|) has an 
asymptotic expansion of the form 

log[Z N (t 3 ,t 4 , . . . ,t 2v )/Z N (0)] = iV 2 e (< 3 ,£4, • ■ ■ Mu) + ei(t 3 ,i 4 , . . . ,t 2u ) + ... 

where the coefficients e s (t3, t±, . . . , t2v) are generating functions for counting labeled, connected, oriented 
maps embedded on a Riemann surface of genus g. 

A map D on a compact, oriented connected surface X is a pair D = (K(D), [t]) where 

(1) K(D) is a connected 1-complex; 

(2) [i] is an isotopical class of inclusions b : K(D) — ► X: 
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(3) the complement of K(D) in X is a disjoint union of open cells (faces); 

(4) the complement of Kq{D) (vertices) in K{D) is a disjoint union of open segments (edges). 

The number of edges connected to a vertex z 6 Ko(D) is called the degree of z. We may represent maps as 
fat or ribbon graphs embedded on an oriented connected surface X. See also the theory of dessins d'enfants 

EDI- 



Ercolani and McLaughlin [5] rigorously showed that log(Zjv) 
form 



possesses an asymptotic expansion of the 



ZN(t3,t4, ■ ■ ■ ,hv) 



1 



— eo(t3)*4) ■ ■ ■ j hu) + -j^2 e i(h, t 



tuv) + 



. , t 2v ) are analytic in a neighborhood of 



Z N (Q) 

in a non-trivial (t^, t^, . . . , tiv) domain. The functions e. g (t%, t^, . 
and are counting functions for labeled maps. 
The function 

oo 

e g (t 3 ,t4, ■ ■ ■ ,t 2v ) = ^2 K g(J^h, 

where n g {j3,jA, ■ ■ ■ ,j2v) is the number of labeled, genus g, maps with Kq(D) containing: j'3 vertices of degree 
3, j4 vertices of degree 4, . . . , j 2u vertices of degree 2u . 

Harer and Zagier [5] solved this enumeration problem in the special case of 1-vertex maps (or monopoles). 
Let G n (N) be the generating function for oriented maps with a single vertex, that is 



f 32v 



G n (N) = N 



1+n £< 

g>0 



where o„ i9 is the number of oriented maps with 1 vertex, and n edges, which are embedded into a genus g 
surface. This generating function is a finite series in iV and is given by the matrix integral 

1. 



(1.3) 



G n {N) = A N / Tr(M^)exp 



H, 



-Tr(M ) 



dAL 



where the integral is taken over the space of N x N Hermitian matrices, and 

-1 



A N = 



exp 



1 



Tr(M 2 ) 



dM 



Theorem 1.1 (Harer-Zagier 8 ). Harer and Zagier explicitly evaluated flJip . to find that 



G n (N) 



(2n)\ 



2™n! ^ 

k=0 



N 
fc + 1 



Bauer and Itzykson show that the data of a labeled map is equivalent to a pair of permutations 
satisfying some compatibility conditions. This result gives an efficient method to determine the genus of a 
map given by its pair of permutations. What emerges is an algorithm for counting the number of genus 
g maps with a particular structure on Kq(D). The idea is to count the number of pairs of permutations 
satisfying the characterization which produce a genus g map. Zvonkin 13 j expanded on this idea. 

The particular partition function Zjv(0,0, . . . , tzv) has been the subject of an extensive calculation UJ. 
The conclusion of this calculation is a prescription for calculating e a (0, 0, . . . , t<iv). Explicit formulas for these 
functions are found for low values of g (g = 0, 1, 2, 3) in terms of an auxiliary function. In the case of g = 
this result produces a closed form expression for the Taylor coefficients Kq 2 ^ (0, 0, ... , j2u)- In the other cases 
worked out explicitly this result is used to find the first few Taylor coefficients of each of the functions. 

A purely combinatoric method for producing an expression for eo(i3,i4, . . . , t 2v ) in terms of generating 
functions of fundamental objects is given in 

The algorithm detailed in this paper has been developed as an independent method for checking the 
detailed work on the fine structure of e g (0,0, . . . , t^J) carried out in and of deriving the constants of 
integration for the formulas derived in that work, and in the author's doctoral thesis. As of the writing of 
this paper the calculation has been done for the Taylor coefficients satisfying vj < 20 and we have found 
exact agreement for all of these terms. 

We will conclude this paper with a modification of our method for counting the number of unoriented 
maps or Mobius graphs. These are maps which are embedded into an unoriented surface. To be more precise 
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the objects we will count are Mobius graphs, ribbon graphs where the edges are allowed to twist. This is 
the setting of interest if the partition function l|1.2fl is replaced with 

'1 



(1.4) 



Z^\t 3 ,t4, . . . ,t 2u ) — 



/ exp 


-NTr ^ 


>s N 





+ t 3 M s + t 4 M 4 + ■■■+ t 2 „M 



dM, 



where the integral is over N x N symmetric matrices. This function is called the partition function of the 
Gaussian orthogonal ensemble (GOE). 



In this case the hypothesis is that 



%N (*3,i4) ■ ■ ■ 7 hu) 



Zn(0) 

where (at least formally) 



Eo(t3, t±, . . . , t 2v ) + —Ei(ts,t4, 



,t 2 u) + 



, t 2l 



u 



K x (j 3 ,j 4 ,...,j 2 »)t3t{ 



■ L 2v ' 



Euler characteristic unoriented maps with Kq(D) 
j2v vertices of degree 2v. The justification of 



i t 2 v) — 

where K x {ja, j±, . . . , j 2v ) is the number of labeled, 
containing: j$ vertices of degree 3, ji vertices of degree 4, 
this enumeration problem is given in □ and 0. 

The existence of the expansion (|1.5I) is not rigorously shown to date. In particular a proof is needed of 
the analyticity of E x . We will give a heuristic argument for the counting property of E x . 

Goulden and Jackson [7] generalized the calculation of Harer and Zagier |S] for 1-vertex maps (monopoles). 
Let F n (N) be the generating function for unoriented maps with a single vertex, that is 



F n (N) = N n ' 



X<2 



where f n<x is the number of unoriented maps with 1 vertex, and n edges, which are embedded into a genus 
g surface. This generating function is a finite series in TV and is given by the matrix integral 

1. 



(1.6) 



F n (N) = B 



N 



Tr(ikP n )exp 



: Tr(M 



dM, 



where the integral is taken over the space of N x N symmetric matrices, and 

-i -l 



Bn — 



exp 



-jTr(M 2 ) 



dM 



Theorem 1.2 (Goulden-Jackson 0). Goulden and Jackson explicitly evaluated i)i.6|) . to find that 



A:=0 



r=0 



2 n n\ ^ 

k=0 



N - 1 
fc + 1 



Our algorithm agrees with the numbers given by F n (N) for all cases checked. What we do not have is a 
method for generalizing Theorem ll.2l to larger families of unoriented maps. 

In section|2]we outline the heuristic argument which motivates using the terms of the asymptotic expansion 
of the partition function H1.2|l as the generating functions of the number of maps of genus g. In section wo 
illustrate the result of pQ which gives a one-to-one correspondence between maps and pairs of permutations, 
of a restricted type. So the group structure provides an algorithm for counting the number of maps of a 
specified type. Section 01 is a collection of tables summarizing some of the calculations we have done with 
this algorithm. It is provided as a means of checking with other known results such as the 1-face maps jS] , 
[T2"| and those in [2] and [Hj. In addition the tables contain all the numbers which were needed to evaluate 
the constants of integration found in . In section we redo the heuristic argument of section in the 
setting of an integral over the symmetric matrices rather than the Hermitian matrices. We show that the 
equivalent counting problem we are concerned with is that of counting maps embedded into an unoriented 
surface. In section [S] we generalize the results of section [3] and associate to each unoriented map a triple of 
permutations. This gives a one-to-one correspondence between unoriented maps and triples of permutations 
of the given type. 
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2. The Wick lemma and Hermitian matrix integrals 



We will first outline the heuristic argument that the e g (ts,t4,, . . . ,t2 U ) are counting functions of labeled 
genus g maps (see [2] and ^2]). This calculation centers around the Wick lemma for Gaussian expectation's: 
The Wick lemma states 

Lemma 2.1. If (•} is a Gaussian expectation, and l\, 1%, . . . , li n are homogeneous linear functions then 

(hh ■ ■ ■ hn) = y^ijji) (^2^2) • • ■ (lijj n ), 
where the sum is taken over the set 1 = i\ < 12 < • • ■ < i n < 2n and ik < jk- 
We are working with the Gaussian Expectation 



(f(M)) = C N /(M)exp 



-N-Tr (M 2 ) 



dM, 



where the integral is taken over the space of N x N Hermitian matrices, and Cjv is a normalizing constant. 
The quadratic expectations of this probability are 



(2.1) 

Consider the partition function 
(2.2) Z N {U) - C N 



/ exp 


-NTi ^ 







l/f2 



1 r4 



dM. 



We may expand out the t^ expression as a Taylor Series and (formally at least) commute this sum with the 
integration to write 



(2.3) 



00 ( 4. \n 

Z N (U) = J2 [ - J £ L N n {[ r n(M 4 )] n ). 



71=0 



Write H2.3fl as a function of the entries of M: 



(2.4) 



( II M imj m M J m km Mk m l m M lmim ) 



1 + i-U) (MijM jk M kl M H ) + 



i,j,k,l 

If we use Wick's lemma on the n = 1 term of 1|2.4[) we find that the contribution is given by 

(2.5) {M l3 M jk ){M k iM u ) + (MijM k i)(Mj k Mu) + {M l3 M u ){M jk M k i) 

Each of these terms is non-zero if i,j,k, and I satisfy the conditions for the corresponding quadratic 
expectations to be non-zero. In this example we find the conditions are respectively i = k, i = I = k = j, 
and j = I. In an analogy to Feynman diagrams these reductions are cataloged by families of maps (see for 
example Figure^. 



Moreover the number of relations between the indices is the number of faces of the corresponding map, while 
the contribution to the expectation (exp [— iVTr (-M 4 )] } is given by the Euler Characteristic of the map. 

The fact that the maps organizing the contributions to the Hermitian matrix integral 12.21) are oriented 
is because of l|2.1f) . 

The final step is to take a logarithm of Z^(t^) to restrict our consideration to connected maps. This 
argument easily generalizes to include other times. 
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Figure 1. The three ribbon graphs representing the maps associated to the n = 1 term of 1)2. 4[l . 



3. Group data associated to each map [D] 

We will now outline the group data associated to a map, the data comes in the form of a pair of permu- 
tations. These two permutations provide a mechanism for counting the number of maps by computing the 
number of permutations satisfying some conditions. Conveniently this permutation structure also provides 
an efficient method of computing the structure of the faces of the map and therefore of calculating the Eulcr 
Characteristic of the map. 

Let E be the number of edges of K(D). We will view the permutations of S2E as acting on the set of 
darts of K{D), which is represented by 

Sl = Sl!U Q 2 , 

where 

Qi = {(u, e) : v is a vertex and e is an edge with two distinct vertices one of them being v } , 

and 

^2 = {(«, e, ±) : v is a vertex and e is an edge with a single vertex v} . 
The element (v,e) G f2i represents the dart based at v and going along e. The element (v, e, ±) G f^2 
represents the dart based at v going along e in the counterclockwise (resp. clockwise) orientation. For each 
edge there are two darts therefore fi = 2E and we can think of S2E as acting by permutations on the set 

n. 

Given a map, (K(D), [t]), we define a subgroup of S2E generated by two permutations (a, r). The 
orientation on X induces (via [t]) a cyclic ordering on the darts attached to each vertex; the first permutation 
a is given by this action. Explicitly a maps the clement (v, e) G ill to the clement (v, e) G f2i or (v, e, ±) G 
where e is the edge counter clockwise in the orientation at v from e. Likewise, a maps the element 
(v,e, ±) G f^2 to the element (i>,e) G fii or (v, e, ±) G where e is the edge counter clockwise in the 
orientation at v from e. The second permutation t, is given explicitly as the permutation which acts on 
Hi by sending (v, e) to (5,e) where 5 is the other endpoint of e; and r acts on Q.2 by sending (t>,e,±) to 
(v,e, =f). 

One sees that c is a product of independent cycles; the length of the cycles in a correspond to the degrees 
of the vertices in Kq(D). The permutation r is a product of £7 disjoint 2-cycles. 

The orbit of < a, r > (the group generated by cr and r) of any element in fl, is all of 17 because K(D) is 
connected. 

To compute the genus of a map from (a, r) we compute the genus of the associated graph from the Euler 
Characteristic 

(3.1) x^V-E + F; 

where V, the number of vertices, is the number of cycles in a; where E, the number of edges, is 1/2 the 
number of elements acted upon by a and r, or is the number of 2-cycles in r; and where F, the number of 
faces, is the number of cycles of the permutation a or (this is merely computing each face by finding all the 
edges which border it). 
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This leads to the theorem: 



Theorem 3.1 (Bauer and Itzykson pQ). There is a one-to-one correspondence between connected maps 
D = (K(D), [l]) and pairs of permutations (a, t), where a is a product of disjoint cycles, and t a product 
of E disjoint 2-cycles, which satisfy connectedness. 

To complete the proof we show that given (a, r) as above one may determine D. The condition of 
connectedness is that the orbit of the group (a, t) ■ 1 is all 2E letters. The following steps produce a map: 

(1) Let K (D) be the set of cycles of a; 

(2) Let K\(D) be the set of 2-cycles of t; 

(3) Each 2-cycle in r ties together two points of Kq(D) identifying the two endpoints of that element of 
Ki(D): 

(4) The Euler Characteristic gives the genus of X ; 

(5) Each cycle in a induces the orientation on X in a neighborhood of l(Kq(D)), this is sufficient to 
determine both X and i up to isotopy (see section 0. 

The punchline is that what we now have is cin algorithm for computing fcgi^JZi J*4 ? ■ ■ ■ )j2v 

) for finite g, j 3 , 

j 4 , • ■ ■ , hu- Let 

2E = 3j 3 +4j 4 + --- + 2vj 2v . 

Fix a to be a permutation formed by j 3 cycles of length 3, j± cycles of length 4, . . . , j 2v cycles of length 2v 
formed from a permutation of n letters. Then we choose each product of E disjoint 2-cycles in Sze- Check if 
(a, t) is connected (by verifying that the orbit of < er, r > -1 is all 2E letters). If (a, r) is connected compute 
the genus as above H3.1[) with \ = 2 — 2g and bin the result. 

3.1. Example . As an example we will show how this algorithm works by computing K g (j 3 — 2,ji = 1). 
Fix a = (1 2 3) (4 5 6) (7 8 9 10). Then choose each r which is a product of 2-cycles of {1,2,..., 10}. 

For example, we might choose r = (1 2)(3 4) (5 6)(7 8)(9 10); for this choice of r the orbit of < a, t > 
forms two disjoint sets {1,2,..., 6} and {7, 8, 9, 10}. 

Choose r = (1 2) (3 4) (5 8) (6 7)(9 10); for this choice of r the orbit of < cr,r > is the set {1,2,. .. , 10}, 
therefore (<r, r) represent a map. We compute the genus of this map by computing the number of distinct 
orbits of a-r. We find that a ■ r = (1 3 5 9 7 4) (2) (6 8) (10), therefore the Euler characteristic is 3-5 + 4 = 2 
and the genus is 0. 

Choose t = (1 2) (3 4) (5 7) (6 9)(8 10); for this choice of r the orbit of < ct,t > is the set {1,2,. .. , 10}, 
therefore (a, r) represent a map. We find that a-r — (1 3 5 8 7 6 10 9 4) (2), therefore the Euler characteristic 
is 3 — 5 + 2 = and the genus is 1. 

4. Results 

In tables HP we present a selection of the results of this calculation, included are the numbers which were 
needed for the calculations in UJ. 

The algorithm can also be used to compute the number of maps of mixed type. For the example from 
section |3~T1 of 2 vertices of degree 3 and one of degree 4 we find 432 genus maps and 468 genus 1 maps. 
For the number of maps with one vertex of degree 3, one of degree 4, and one of degree 5: 2160 (genus 0), 
6480 (genus 1), and 1440 (genus 2). 

5. Wick lemma and symmetric matrix integrals 
This section is motivated by the partition function given by \1A\ . 

We will now show, using the Wick lemma, that the asymptotic expansion of log ^zj^^ enumerates 

unoriented maps. This will be a formal calculation only, the necessary analyticity of the terms of this 
expansion has not been rigorously shown. We present the calculation here as a motivation for the counting 
problem we are computing. 

We are working with the Gaussian Expectation 

(f(M)) = C N [ /(M)exp r 



-N-Tt(M 2 ) 
4 



dM, 



Table 1. First the number of one vertex maps found by [Hj (whose duals are the well 
studied 1-face maps |12| ~1 



Degree 


Genus 





1 


2 


3 


4 


5 


4 




2 


1 










6 




5 


10 










8 




14 


70 


21 








10 




42 


420 


483 








12 




132 


2310 


6468 


1485 






14 




429 


12012 


66066 


56628 






16 




1430 


60060 


570570 


1169740 


225225 




18 




4862 


291720 


4390386 


17454580 


12317877 




20 




16796 


1385670 


31039008 


211083730 


351683046 


59520825 



Table 2 . The number of two vertex maps 



Degree 


Genus 





1 


2 


3 


4 


3 




12 


3 








4 




36 


60 








5 




180 


600 


165 






6 




600 


4800 


4770 






7 




2800 


34300 


81340 


16695 




8 




9800 


215600 


1009400 


781200 




9 




44100 


1323000 


10478160 


19158300 


3455865 


10 




158760 


7408800 


94091760 


333774000 


218402730 



Table 3. For the case of degree 4 maps we sum up our calculations in the following table 
(these numbers agree with those found in [5] and 0) 



Vertices 


Genus 





1 


2 


3 


1 




2 


1 






2 




36 


60 






3 




1728 


6336 


1440 




4 




145152 


964224 


770688 




5 




17915904 


192098304 


348033024 


58060800 



where the integral is taken over the space of N x N symmetric matrices, and Cm is a normalizing constant. 
The quadratic expectations of this probability are 



(5.1) 

Consider the partition function 



{MijMki) = — [SuSjk + SikSji] 



/ exp 


-NTi ^ 


>s N 





dM. 



We may expand out the expression as a Taylor Series and (formally at least) commute this sum with the 
integration to write 



Write this expression as a function of the entries of M: 

Z { N ] {t A ) = ^^-N n Yj ( II M ^J m M lmkm M kmlm M lml J 

= 1 + (-*4) (Mi,jM j>k M ki iMi ti ) + ... 

i,j, k ,l 

If we use Wick's Lemma on the n = 1 term we find that the contribution is given by 

(5.2) (M l3 M jk ){M k iM H ) + (MijM k i){MjkMu) + {M l3 M kl )(M, k M k i) 

Each of these terms is non-zero if i, j, k and I satisfy the conditions for the corresponding quadratic expec- 
tations to be non-zero. As in the case of oriented maps these reductions may be cataloged by families of 
maps. In this case the maps are allowed to have edges that reverse orientation because of formula H5.ll) . 
Focus on the middle term of (|5.2|) 

(5.3) {M ij M kl ){M jk M u ). 

There arc four cases which give that this contribution is non-zero: i = l = k = j,i = k = l= j,i = l = k= j, 
or (i = k,j = I), these correspond to the four unoriented maps represented by the Mobius graphs in Figure 

El 




Figure 2. Four Mobius graphs corresponding to the unoriented maps whose underlying 
oriented map is the single, genus one, oriented map with one vertex of degree four. 



The final step is to take a logarithm of Z N {t^) to restrict to connected unoriented maps. 

6. Unoriented maps 

We will lift the algorithm for counting labeled connected oriented maps to count labeled connected un- 
oriented maps. We will first go over the modifications which must be made to the group data to account for 
the independence that the local orientations of each vertex have. 
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An unoriented map on a compact unoriented connected surface X is a pair D = (K(D), [t]) where 

(1) K(D) is a connected 1-complex; 

(2) [l] is an isotopical class of inclusions l : K(D) — ► X; 

(3) the complement of K(D) in X is a disjoint union of open cells (faces); 

(4) the complement of Kq(D) (vertices) in K(D) is a disjoint union of open segments (edges). 

In the same way that ribbon graphs are associated to an oriented map we associate a Mobius graph (ribbon 
graph with orientation reversing and orientation preserving edges) to an unoriented map. 

Given an unoriented map D, there are three permutations associated to the data contained in D. To 
begin, for each vertex, v e K n (D) choose a local orientation, ri v induced by the embedding l. We will view 
the four permutations as acting on the set of darts xZ 2 , given explicitly as: Q — fli U Q2, where 

O x = {(v, e, ±n v ) : v is a vertex and e is an edge with two distinct vertices one of them being v} , 

and 

^2 = {(v, e, ±n Vl ±) : v is a vertex and e is an edge with a single vertex v} . 

The element (v, e, ±n„) £ Hi represents the dart based at v pointing along e together with the local orien- 
tation ri v (resp. — n v ). The element (v,e, n v ;±) E represents the dart based at v which traverses e in 
the counter clockwise direction with respect to the local orientation (resp. clockwise direction with respect 
to the local orientation) together with the local orientation ft v . Likewise the element (v,e, — n v ;±) G 2 
represents the dart based at v which traverses e in the counter clockwise direction with respect to the local 
orientation (resp. clockwise direction with respect to the local orientation) together with the normal vector 
-n v . 

An edge is said to be orientation reversing if n v at v on the edge is — nj, after translating to v along the 
edge; otherwise we say that the edge is orientation preserving. This is a property of both the map and the 
choice of local orientations. 

We find three permutations (0, a, t) : 

On fli we define the action of (cf>, ct, t) by: 

• (p acts by sending (v, e, ±n v ) to (v, e, ^n v ). 

• a acts by sending (w, e, n v ) (resp. (w, e, — n v )) to the dart which is counter clockwise (resp. clockwise) 
from (v,e) at v, maintaining the orientation ri v (resp. —ft v ). 

• t acts by sending (v, e, ±n„) to (v, e, ±n v ) when e is orientation preserving or (v, e, ^fn v ) when e is 
orientation reversing. 

On f^2 we define the action of (</>, a, r) by: 

• <p ac ts by sending (v, e, ±n v ; ±) to (v, e, ^n v ; ±). 

• g acts by sending (v,e,n v ;±) (resp. (v, e, — ri v ; ±)) to the dart which is counter clockwise (resp. 
clockwise) from (v,e, ±) at v, maintaining the orientation ri v (resp. — n v ). 

• t acts by sending (v,e,n„;±) to (w, e, ±n„; =p) (resp. sending (w, e, — ri v ; ±) to (w, e, =Fn„; =f) ) de- 
pending on whether r is orientation preserving or reversing 

From this definition the following properties are apparent: 

(1) (f)(T(f) = (T _1 and <pT<p = t 

(2) (j) an d t are fixed point free permutations, each is a disjoint product of 2-cycles. No 2-cycle of r is a 
2-cycle of 4>. 

(3) a is a disjoint product of 2j3 3-cycles, 2j4 4-cycles, . . . , 2j2 V 2^-cycles. 

(4) A connected map has the property that (0, a, r) has a single orbit. 

(5) t contains the data of how the darts of K(D) are attached together to form the map; together with 
<p, t determines how the local orientations relate globally. 

(6) The cycles of a o t correspond to faces of the map. There are two cycles of equal length for each 
face, representing the two directions the face can be traced along its edges. 

This leads to the theorem: 

Theorem 6.1. There is a one-to-one correspondence between unoriented maps and triples {<p, a, t) satisfying 
conditions (1-6). 
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We will complete the proof by showing that given three permutations ((f), a, r) satisfying conditions (1-6), 
we can reconstruct (up to the choice of local orientation) the map (K(D), [4]) and the surface X. 

A. Let K (D) be the set of cycle pairs (0-1,0-2) in a such that <$>a\§ — a 2 x . 

B. From each pair (oi, <r 2 ) choose <j\ to represent the counter clockwise or up orientation. 

C. Let Ki(D) be the set of 2-cycle pairs (t 1 ,t 2 ) in r such that 4>Tx<j) = t 2 . Each pair (t 1 ,t 2 ) connects two 
(possibly the same) elements of Kq(D). If (ti,t 2 ) connects the up orientation to the up orientation this 
edge is orientation preserving, otherwise it is orientation reversing. 

D. To construct the surface X and the embedding class [4] do the following: 

(i) For each vertex (0-1,0-2) draw a polygon with sides labeled and oriented by the elements that o\ 
acts on. 

(ii) For an orientation preserving edge connecting two vertices, glue together the corresponding sides of 
their polygons so that the relative orientations agree at the corresponding side. For an orientation 
reversing edge connecting two vertices, glue together the corresponding sides of their polygons so 
that the relative orientations are opposite at the side. 

In conclusion: up to the choice of local orientations at each vertex, the data given by (</>, a, r) produces 
a map. Likewise, a map, up to the choice of local orientations at each vertex, produces a triple ((f), a, t) 
satisfying conditions 1-6. 

6.1. The algorithm. We now have all the tools necessary to lift our algorithm to the enumeration of 
unoriented maps. 

Fix to be a fixed point free disjoint product of 2-cycles. Fix a to be a product of 2j3 3-cycles, 2j4 
4-cycles, . . . , and 2j2!/ 2^-cycles, such that (f> connects cycles in a into pairs (oi, 0-2) of the same cycle type, 
choose o-i to represent the ^( CT i,ct 2 ) orientation. 

Cycle through all fixed point free products of disjoint 2-cycles r such that conditions 1-6 are satisfied 
by (0, o-, t), in particular no 2-cycle in (f> appears in r. If, for this choice of r, the corresponding map is 
connected, then compute the Euler characteristic: the number of vertices is half the number of cycles in a, 
the number of edges is half the number of 2-cycles in t, and the number of faces is half the number of cycles 
in a ■ t. Bin the results sorted by the Euler characteristic. 

7. Conclusion 

The algorithm presented here for counting oriented maps was used in [^] and the authors dissertation 
to evaluate constants of integration appearing in the explicit calculation of e g (t 2v ) for low genus (g < 3). 
Additionally the algorithm provided a useful method for verifying the validity of the results in 

In looking beyond the calculations carried out in |S] it became apparent that it would be useful to have 
a similar algorithm to compute the number of unoriented maps. The permutation structure associated to 
oriented maps has a precise analogue in unoriented maps. This method may prove useful in further studies 
of the GOE and GSE matrix integrals, and associated problems. 
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